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Abstract 

We use molecular simulations to explore how sample dimensions and interfacial 
properties impact some generic aspects of the mechanical and structural behavior 
of nanoconfined materials. Specifically, we calculate the strain-dependent proper- 
ties of minimum-energy thin-film particle configurations (i.e., inherent structures) 
confined between attractive, parallel substrates. We examine how the relationship 
between the transverse strain and the stress tensor (the equation of state of the 
energy landscape) depends on the properties of the film and substrate. We find that 
both film thickness and film-substrate attractions influence not only the mechanical 
properties of thin films, but also the shape and location of the "weak spots" where 
voids preferentially form in a film as it is strained beyond its point of maximum ten- 
sile stress. The sensitivity of weak spots to film properties suggests that nanoscale 
materials may be intrinsically vulnerabile to specific mechanisms of mechanical fail- 
ure. 
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1 Introduction 



Materials confined to very small spatial dimensions behave differently than in 
the bulk. In addition to showing quantum confinement effects, they also dis- 
play thermodynamic, kinetic, and mechanical limits of stability that depend 
on sample size, shape, and the characteristics of their interfaces. Specific exam- 
ples of property modifications include surface-induced phase transitions, shifts 
of t he bulk glass transition, and interface-mediated modes of mechanical fail- 



ure (Gelb et. al . 



1999; 



Forrest and Dalnoki-Veress 



2001 



Hutchinson and Suo 



19921) . Unfortunately, because molecular-scale processes in highly inhomoge- 
neous environments are difficult to resolve experimentally, a mechanistic pic- 
ture for precisely how nanoconfinement impacts stability has been slow to 
develop. This presents a practical barrier to the design of technological ap- 
plications, in particular those relying on solid-state nanostructures to exhibit 
mechanical integrity over a broad range of conditions. 

In this Article, we study an elementary model system that sheds new light on 
how sample dimensions and interfacial properties can influence the mechanical 
behavior of nanoconfined materials. Specifically, we use molecular simulations 
to calculate the strain-dependent properties of mechanically-stable films of 
particles confined between attractive, parallel substrates. Although analogous 
studies have been carried out for models of isotropic materials, this is, to 
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our knowledge, the first systematic investigation of the relationship between 
the transverse strain and the pressure tensor of the inherent structures (min- 
imum potential energy configurations) of highly inhomogeneous films. Our 
main finding is that both sample dimensions and substrate attractions sub- 
stantially influence not only the mechanical properties of thin films, but also 
the morphology and location of "weak spots" where voids preferentially form 
in a film as it is strained beyond its point of maximum tensile stress. Although 
the precise role that these weak spots play in dynamic deformation processes 
is presently unknown, they appear intimately linked to material failure by 
quasistatic tensile deformation. Moreover, the sensitivity of weak spots to film 
properties suggests that nanoscale materials may be intrinsically vulnerabile 
to specific mechanisms of mechanical failure. 



Since plastic deformation and failure are inherently dynamic events, and since 
the molecular-scale rearrangements that accompany them in amorphous mate- 
rials are still poorly understood, molecular dynamics (MD) simulations would 
appear to represent an ideal theoretical tool for their investigation. In fact, 
MD simulation studies over the past decade have been instrumental in gaining 



(Falk and Laneer 




1998 
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GersaDDe and Rc 


tbbins 




1999 


; E 


cottier and Robbins. 


2001 




200a 


StevenJ. 


2001 


GersaDDel. 


2002 


CaDaldi et al. 


200a 


Varnik et al. 


2004 


Van Workum and de Pablo 


■ 


200,4 


^oshimoto et a 


L, 


200^ 


). These in- 



sights have facilitated the interpretation of experime nts and have aided in th e 



1998). 



introduction of simple theories for viscoplastic flow (jFalk and Langer 
On the other hand, despite recent advancements in algorithms for long-time 
dynamics, MD simulations are still limited to accessing relatively short time 



and length scales. Thus, it is currently computationally prohibitive to use MD 
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to exhaustively explore the effects that sample dimensions and interfacial con- 
ditions have on the mechanisms of mechanical failure, even for simple model 
systems. The development of alternative methods for probing the molecular- 
scale origins of failure in glasses is still essential. 

One complimentary approach is to calculate how the properties of a mate- 
rial's mechanically-stable inherent structures depend on the state of strain. 
This strategy has typically been implemented in one of two ways to investi- 



gate defor mation, ten s ile strength, and fai l ure of 



method 



Vlott et al 



1993 



Hutnik et al 



1993 



amorphous solids. T 



re first 



Malandr o and Lacks!. 



1997 



1998 . 



1999; 



Lund and Schuh. 12003a b 



Malonev and Lemaitre 



2004aU rD calls 



for subjecting an ensemble of inherent structures, created at a prescribed state 
of strain or stress, to an athermal quasistatic deformation process consisting 
of alternating steps of small affine deformation follow ed by potential energy 



minimization. 



Malandro and Lacks! ( 1997 



1998 . 



1999) have used this protocol 



to investigate the connection between strain-induced plastic rearrangements 
in amorphous materials and the annihilation of minima on the material's po- 
tential energy landscape. A s imilar implementation has also been used by 



Malonev and Lemaitre 




2004 



to study the energy fluctuations associated 
with amorphous plasticity and the behavior of elastic constants near the onset 
of material failure. 



Alternatively, one can generate collections of inherent structures at each macro- 
scopic strain state of interest by mapping equilibrium particle configurations 
from high-temperature sim ulations to their local potential energy minima 



flStillinger and Weber 



19821 ). This procedure, which we adopt h ere, has been 



prim a rily employed 



19J7J 



) V Deb e nedetti, St i 



Roberts et al 



1999; 



Htz et al. 



linge r , and collaborators (pastry et al 



2001 



Shen et al 



20021 ) to determine 
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how the inherent structure pressure of 'bulk' glass-formers depends on density, 



a relationship t hat has been termec 



scape (EOSEL) ([Debenedetti et al 



the equation of state of an energy land- 



1999). Several trends have emerged from 



simulated EOSELs that give insights into the mechanical properties of amor- 
phous solids. For example, inherent structures formed with densities higher 
than the material's Sastry density p a (shown in Fig la ) are structurally ho- 



mogeneous, whereas "weak spots" (pastry et al 



19971 ) contained in lower- 



density equilibrium configurations develop into fissures or voids upon energy 
minimization. Thus, ps is a material property that represents the minimum 
density for which mechanically-stable solid structures can remain structurally 
homogeneous and void free. Moreover, the corresponding isotropic tension 
—pis{ps) is the maximum amount that an inherent structure of that material 
can sustain prior to rupture. Note that ps obtained from the EOSEL is concep- 
tually similar to the density of maximum isotropic tensile stress o btained from 



a qua sistatic expansion process, and, in fact, recent simulations (jShen et al 



2002) show quantitative agreement between these two densities for a number 



of model materials. 



2 Probing Vulnerabilities to Failure in Nanoscale Films 



The EOSELs of ult rathin films have been recently predicted by an app roximate 



mean-field theory (jMittal et al 



2004; 



Truskett and Ganesan . 



20031 ). but, to 



our knowledge, they have never been determined via molecular simulation. 
One naturally expects far richer behavior in thin films than in bulk materials 
because interfacial interactions render them anisotropic and statistically inho- 
mogeneous. One consequence, illustrated in Fig. lb, is the emergence of two 
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distinct versions of the curve shown in Fig. la, representing the density depen- 
dencies of the transverse (p\s) and normal (pfs) components of the inherent 
structure pressure tensor. This additional dimension of the EOSEL prompts 
several new questions pertaining to the film's intrinsic vulnerabilities to me- 
chanical failure. For example, how do film thickness and substrate attractions 
affect the directionality of mechanical failure (i.e., which component of the 
stress tensor will exhibit a maximum at a smaller value of strain)? Moreover, 
is structural failure (i.e., void formation') in i tiated at the strain of maxim al 



stress as it is in bulk materials (pastry et al 



1997 



Debenedetti et al. 



1999)': 



If so, do film thickness and substrate attractions strongly influence the location 
and the morphology of the corresponding voids? 

To investigate these issues, we calculated the EOSEL of model films com- 
prising N = 864 Lennard- Jones (LJ) particles (truncated and shifted at 



2.5 (jSastrv et al. 



19971 )) confined between parallel substrates that in- 



teract with them through an effective 9-3 LJ potential, 



( \ 2lX 3 



15 V z 



Z 



Here, z is the dis tance between substrat e and particle center, p w a, 



0.988, 



and a w = 1.0962 rtEbner and Saam 



19_7_7J). All quantities reported in this work 



are implicity nondimensionalized by the standard combinations of length and 
energy scales provided by the effective diameter a and the well-depth e of 
the LJ film particles, respectively. Periodic boundary conditions were applied 
in the x and y directions to simulate a system of infinite transverse extent. 
Additional simulations (not shown here) using N = 1024 and N = 2048 
particles were also carried out to verify that the results were insensitive to N. 
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The effects of film thickness and strength of the film-substrate attractions 
were examined by simulating nine films characterized by the permutations of 
the following substrate separations L G {5, 7.5, 10} and film-substrate we ll 



19771) . 



depths e w e {0.2ej w , €f w ,5ef w }, where ej w = 1.2771 (jEbner and Saam 
The inherent structures required for the EOSEL analysis were obtained by 
mapping high-temperature (T = 2.5) equilibrium configurations from a se- 
ries of fluid film densities in the range p = 0.2 — 1.0 to their local poten- 



tial energy minima using LBFGS ([Liu and Nocedal 



1989), a limited-memory 



quasi-Newton multidimensional optimization routine. The equilibrium config- 
urations were generated by Monte Carlo simulations in the canonical ensemble. 
For every film type and strain state, 100 equilibrium configurations (separated 
by 500 Monte Carlo cycles) were selected for minimization. The components 



of the inherent structure pressure 



cal 



censo r were calculated both by the global 



viri al equation (IBorn a nd Huang, 



1954 ) and by spatially averaging the lo- 



Irving and Kirkwooc 



expressions. The two approaches produced 



indistinguishable results. 



For the structural analysis of the inherent structure configurations, void space 
was defined to be the volume of film comprising all points that lie both a dis- 
tance d > 1 from any particle center and d > (a w + 1)/2 from either substrate. 
The void probability is then simply the fractional film volume available for 
insertion of an additional particle center without creating "overlap" with an 
existing particle or substrate. For each strain state, the average void probabil- 
ity was calculated by performing 10 6 trial insertions of hard-core test particles 
of unit diameter. These data were binned according to normal position z 
to determine the film's inhomogeneous void probability profiles. The meth- 
ods described above for both generating and analyzing our thin-film struc- 
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tures follow directly from established techniques for studying bulk materials, 
an d the interested reader can find more information in recent review articles 
bv bebenedetti et all jlQQflL boOlh . 



As a final point, we discuss the main rationale for choosing the LJ model for 
our initial investigation. First, it is a microscopic model that shows qualita- 



tively realistic mechanical and structural behavior (pastry et al 



19971 ). but it 



also is simple enough to allow a systematic investigation of the effects of inter- 



faces and confinement on its p roperties. Moreover , since the bulk L 



is already well characterized (pastry et al 



1997 



Debenedetti et al 



system 



1999) 



any confinement-induced deviations should be relatively straightforward to 
recognize and interpret. Finally, a knowledge of these deviations will provide 
a useful basis for understanding future simulations on materials with richer 
molecular attributes. 



3 Results and Discussion 



In this section, we study how the specific properties of our model LJ films im- 
pact their vulnerability to various modes of mechanical failure. To accomplish 
this, we investigate how sample dimensions and film-substrate attractions af- 
fect the strain dependencies of the stress tensor and the void-space morphology 
of thin-film inherent structures. The locations and shapes of the voids or 'weak 
spots' that appear in the films when they are strained beyond their maximal 
stress states provide markers for the corresponding mechanisms of failure. We 
discuss how these mechanisms arise from structural inhomogeneities in the 
film and the properties of the film-substrate interfaces. 



S 



Fig. 2a shows the simulated EOSEL of an ultrathin film (L = 5) confined 
between 'weakly attractive' substrates (e^ = 0.2ef w ). Note that pj~ s exhibits a 
minimum at a higher value of density than does p\ s , or, equivalently, the max- 
imum normal stress corresponds to a smaller strain than does the maximum 
transverse stress. Thus, from a mechanical perspective, the inherent structure 
film can be considered more vulnerable to failure in the normal direction when 
subje cted to plane strain. Fig. 2a also illus t rates that, as is observed in the 



bulk (pastry et al 



1997 



Debenedetti et al. 



1999), the probability of finding 



voids in the film becomes non-neglible for densities below the maximal stress 
state. This connection between mechanical and structural failure is examined 
more closely in the void probability profiles and film images of Fig. 2b and 
2c, respectively. As can be seen, the main structural effect of transverse strain 
is the formation of planar voids where the film detaches from one of the sub- 
strates, indicating an adhesive mechanism for mechanical failure. 



One natural question to ask is, do substrate attractions noticeably affect the 
film's mechanical and structural vulnerabilities to failure? Fig. 3 shows the 
EOSEL for an ultrathin film (L = 5) with 'strongly attractive' substrates 
(e w = 5ef w ). Comparison of Fig. 2a and 3a shows that, at least for ultrathin 
films (L = 5), there seem to be some features of the EOSEL that are indepen- 
dent of the strength of the substrate attractions. For example, the film confined 
between strongly attractive substrates also shows mechanical vulnerability to 
failure in the normal direction, with pf s exhibiting its minimum at a higher 
value of density than p\ s . Moreover, Fig. 3a shows a similar correspondence 
between the attainment of maximal stress and the initial appearance of voids 
in the film. However, the key difference can be found in the location of the 
'weak spots' where strain-induced voids appear. As is evidenced by the void 
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profiles and film images of Fig. 3b and 3c, the strongly attractive substrates 
promote a cohesive rather than adhesive failure mechanism, indicated by the 
formation of planar voids between the molecular layers in the center of the 
film. This mode of failure results in films with much larger tensile strength 
and toughness than films that fail adhesively. 

The main physical result from Fig. 2 and 3 is that strain-induced local stresses 
in highly inhomogeneous films, dictated by the relative strengths of the inter- 
particle and film-substrate attractions, introduce specific types of vulnerabil- 
ities to failure (e.g., adhesive or cohesive failure). These vulnerabilities have 
important implications for the material properties of the films. Interestingly, 
there is also preliminary evidence to suggest that such energetic vulnerabili- 
ties may strongly influence the actual mechanisms of dynamic failure events. 



In particular, re cent MP simulations of 



j olvm eric thin films confined between 



solid substrates ([Gersappe and Robbing 



1999) have found that the transition 



from adhesive to cohesive mechanisms of failure is essentially independent of 
temperature and, moreover, that it is primarily controlled by the strength of 
substrate attractions. The fact that both dynamic and EOSEL-type simula- 
tion studies provide a similar physical picture is intriguing. However, more sys- 
tematic investigations of the energetic and entropic contributions to dynamic 
failure events in model systems will be necessary to assess the generality of 
these observations. 

To explore the impact of thickness L on the mechanical and structural prop- 
erties of nanoscale films, we consider two additional cases. The first, shown in 
Fig. 4a, is the EOSEL of a thicker film (L = 10) confined between strongly 
attractive substrates (e w = 5ej w ). Note that considerable asymmetry between 
pj~ s and p\ s is still apparent in thicker films, with the former component at- 
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taining a larger magnitude of maximum tensile stress at a slightly smaller 
value of plane strain. This behavior reflects the fact that film particles are 
strongly attracted to the substrates, making it favorable for the film to acco- 
modate plane strain by tearing itself apart internally. However, the proximity 
of the strains of maximal normal and tensile stress in Fig. 4a (as compared 
to Fig. 3a) suggests that one can expect differences between the structural 
mechanisms of failure for thicker versus ultrathin films. 



To interrogate the microscopic origins of these differences, we look to the void 
space analysis. As can be seen by comparing Fig. 4b and 4c to Fig. 3b and 3c, 
cohesive failure does show qualitatively different structural consequences in 
ultrathin and thicker films. While ultrathin films fail by forming planar voids 
between molecular layers parallel to the substrates, thicker films generally 
fail by forming voids that are 'local' and more isotropic (more balanced in the 



transverse and normal directions), not so differ ent from the 1 



jubble-like" voids 



that form from weak spots in bulk materials (jSastrv et al 



1997 



Shen et al 



20021 ) . Note that this local cohesive failure mechanism produces a film with a 



much smaller tensile strength than the ultrathin film of Fig. 3. 



The second thicker-film case (L = 10) that we examine is the one confined 
between neutral substrates (e w = ef w ). What type of failure mechanism is 
expected to prevail for this film? Fig. 4 shows that in thicker films, like in 
bulk materials, there is a strong tendency to cavitate upon expansion due local 
cohesive stresses. However, since the substrates are not nearly as attractive in 
the neutral wall film, it is not clear at the outset whether cavitation will be 
controlled by the cohesive stresses in the normal direction (as in Fig. 4) or in 
the transverse direction. 
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The EOSEL shown in Fig. 5a displays only slight asymmetry in the normal and 
transverse pressure tensor components, indicating that the two aforementioned 
mechanisms for failure (normal cohesive and transverse cohesive) are in close 
competition. In contrast to the results of Fig. 4a, p\ s in the neutral wall film 
exhibits its minimum at a density slightly larger than does pf s . This suggests 
that there is a preference for the transverse cohesive mechanism for failure. 



To see the structural consequences of this type of failure, we examine the 
void probability profiles and the film images of Fig. 5b and 5c, respectively. 
Fig. 5b shows that while voids can appear anywhere in the film, they show 
a modestly higher tendency for forming near the neutral substrates. The film 
configurations of Fig. 5c confirm that the strain-induced voids that form in 
the thicker film are also "bubble-like" like in Fig. 4c, as opposed to the planar 
voids observed in the ultrathin films of Fig. 2 and 3. 



Finally, we point out there is an additional structural consequence of trans- 
verse cohesive failure. Since the initiating voids have a modest preference for 
forming near substrates with weak attractions, and since they grow in size 
upon the application of strain, they can ultimately lead to a secondary adhe- 
sive failure event. This simply means that it may not be easy to distinguish 
between adhesive and transverse cohesive modes of failure "after the fact", 
i.e., by only examining the final structure of the failed sample. In contrast, we 
do not find a similar strain-induced cascading of failure events for thin films 
confined between strongly attractive substrates. 
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4 Summary and Conclusions 



To complete our nanoscale film analysis, the main features of the EOSELs 
for all 9 films studied are summarized in Table 1. The information provided 
includes the film thickness L, the well-depth of the film-substrate potential 
e w/tfw, the direction (|| or _L) of vulnerability to mechanical failure (i.e., which 
stress component shows a maximum at a higher density), and the correspond 



ing density where this o ccurs ps (i.e., the Sastry density (pastry et al 



Debenedetti et al. 



1997 



19991 ) for the film). 



The main trends make good physical sense. For thicker films, we see a crossover 
in vulnerabilities from normal cohesive failure (when confined between strongly 
attractive substrates) to transverse cohesive failure (when confined between 
weakly attractive substrates). In the former case, strong bonding to the sub- 
strates tears the film apart upon the application of strain, whereas the film 
ruptures under its own strain-induced cohesive forces if the substrates are only 
weakly attractive. These purely energetic trends are in quantitative agreement 



with the 



predictions o 



for films ( Mittal et al 



a recently introduced energy land scape-based theory 



2001; 



Truskett and Ganesan 



20031 ). In the thicker film 



cases discussed above, the weak spots where voids form appear "bubble-like", 
similar to bulk material behavior. For ultrathin films, we find a crossover in 
vulnerabilities from cohesive failure (when confined between strongly attrac- 
tive substrates) to adhesive failure (when confined between weakly attractive 
substrates). Films that fail cohesively have much higher tensile strength and 
toughness. These trends are in line with those s een in recent MP simulations o f 



dynamic failure in polymer thin-film adhesives ((Gersappe and Robbins 



1999). 



In the ultrathin film cases, the weak spots appear planar (parallel to the sub- 
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strates), and their location is strongly influenced by the inhomogeneous struc- 
tural layering of the films. 



In spite of the intriguing results provided by the present analysis of the 
EOSELs for model thin films, there are many open issues that warrant fu- 
ture study Even if we confine ourselves to EOSEL thin-film studies, there 
remains a need to understand how factors such as substrate morphology 
and/or the effects of "free" interfaces impact material stability. Neither of 
these issues are addressed in the present work. The question of how molec- 
ular connectivity and architecture (from small molecules to polymers) can 
impact the mechanical strength, structure, and the possible mechanisms for 
failure is also of great fundamental and practical importance. For instance, 
would molecules of ellipsoidal shape show fundamentally different behavior 
from the sph erically-symmetric molecules studied here, since the former have 



been shown 



Donev et al. 



ter (jTorquato et al 



2000; 



200 41) to pack much more efficiently than the lat- 



Truskett et al 



20001 )? Furthermore, how does the 



present picture change if one considers the EOSEL under plane stress rather 
than plane strain conditions? What about failure occuring under shear de- 
formations? Finally, the role that failure mechanisms obtained via quasistatic 
analyses play in actual dynamic failure processes is still far from resolved, and 
it is a issue that demands further scrutiny. 
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Figure Captions 

(1) Schematic of the relationship between inherent structure (IS) pressure 
and density, the so-called equation of state of an energy landscape (EOSEL) 
(a) In bulk isotropic materials, inherent structures formed at densities 
above the Sastry density p$ are structurally homogeneous. Those formed 
at densities below p$, the state of maximum isotropic tensile stress, con- 
tain voids or fissures (Sastry et al., 1997; Debenedetti et al., 1999). (b) 
Schematic of the EOSEL for a thin film, showing the diagonal compo- 
nents of the inherent structure pressure tensor in the transverse (p\g) an d 
normal (pj~s) directions. 

(2) (a) EOSEL of an ultrathin film (L = 5) confined between weakly attrac- 
tive substrates (e w = 0.2ef w ), depicting normal (□) and transverse (>) 
components of the inherent structure pressure tensor. Also shown is the 
corresponding probability of finding of voids (o) in the film, (b) Void prob- 
ability distribution P(z) in the normal direction for p = 0.2, 0.65, 0.66, 
and 1.0. The arrow indicates increasing density (decreasing transverse 
strain), (c) Visual depiction of inherent structure configurations for den- 
sities p = 0.5 (left) and p = 0.66 (right). These illustrate where planar 
voids form near the substrate, indicating an adhesive failure mechanism. 

(3) (a) EOSEL of an ultrathin film (L = 5) confined between strongly at- 
tractive substrates (e w = 5ef w ), depicting normal (□) and transverse 
(>) components of the inherent structure pressure tensor. Also shown is 
the corresponding probability of finding of voids (o) in the film, (b) Void 
probability distribution P{z) in the normal direction for p = 0.2, 0.3, 0.76, 
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and 1.0. The arrow indicates increasing density (decreasing transverse 
strain), (c) Visual depiction of inherent structure configurations for den- 
sities p = 0.5 (left) and p = 0.76 (right). These illustrate where planar 
voids form in between molecular layers, indicating a cohesive failure mech- 
anism. 

(4) (a) EOSEL of a thicker film {L = 10) confined between strongly at- 
tractive substrates (e„, = 5ej w ), depicting normal (□) and transverse (>) 
components of the inherent structure pressure tensor. Also shown is the 
corresponding probability of finding of voids (o) in the film, (b) Void prob- 
ability distribution P(z) in the normal direction for p = 0.2,0.86,0.87, 
and 1.0. The arrow indicates increasing density (decreasing transverse 
strain), (c) Visual depiction of inherent structure configurations for den- 
sities p = 0.6 (left) and p = 0.82 (right). These illustrate where "bubble- 
like" voids form in the film, consistent with a normal cohesive failure 
mechanism. 

(5) (a) EOSEL of a thicker film (L = 10) confined between neutral substrates 
(e w = tfw), depicting normal (□) and transverse (>) components of the 
inherent structure pressure tensor. Also shown is the corresponding prob- 
ability of finding of voids (o) in the film, (b) Void probability distribution 
P(z) in the normal direction for p = 0.2,0.81,0.82, and 1.0. The arrow 
indicates increasing density (decreasing transverse strain), (c) Visual de- 
piction of inherent structure configurations for densities p = 0.4 (left) and 
p = 0.82 (right). These illustrate where "bubble-like" voids form near the 
substrates, consistent with a transverse cohesive failure mechanism. 
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Table 1 

Summary of the main features of the EOSELs for 9 films: film thickness L, the 
well-depth of the film-substrate potential e w /ej w , the primary direction (|| or _L) of 
vulnerability to mechanical failure, and the corresponding density where this occurs 
ps (i.e., the Sastry density for the film). 
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